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We present a measurement of the semileptonic mixing asymmetry for B° mesons, afi, using 
two independent decay channels: B° — > fi + D~ X, with D~ — > K + tt~-k~ "; and B° -!• with 
D*~ — > D°tt~ , D° — > K + tt~ (and charge conjugate processes). We use a data sample corresponding 
to 10.4 fb _1 of pp collisions at ^/s = 1.96 TeV, collected with the DO experiment at the Fermilab 
Tevatron collider. We extract the charge asymmetries in these two channels as a function of the 
visible proper decay length (VPDL) of the B meson, correct for detector-related asymmetries using 
data-driven methods, and account for dilution from charge-symmetric processes using Monte Carlo 
simulation. The final measurement combines four signal VPDL regions for each channel, yielding 
a^i = [0.68 ± 0.45 (stat.) ± 0.14 (syst.)]%. This is the single most precise measurement of this 
parameter, with uncertainties smaller than the current world average of B factory measurements. 

PACS numbers: 11.30.Er, 12.15.Ff, 14.40.Nd 



I. INTRODUCTION 

Fundamental asymmetries in the interactions of ele- 
mentary particles influence the large-scale behavior of the 
universe. Of particular interest is the process of baryogc- 
nesis, whereby an initially symmetric system of particles 
and antiparticles produced by the Big Bang evolved into 
the observed matter-dominated universe of the present 
day. Current theoretical models, building on the work of 
Sakharov [1], require CP-symmetry violating processes in 
order for baryogenesis to have occurred in the very early 
universe [2-5] . As such, studies of asymmetries in parti- 
cle physics experiments have an influence far beyond the 
scale that they probe directly. 

CP symmetry implies that physical processes are in- 
variant under the combined parity and charge conjuga- 
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tion transformations. The standard model (SM) of par- 
ticle physics is not CP symmetric as it stands, due to a 
complex phase in the quark mixing matrix of the weak 
interaction, which has been measured to be non-zero [6] . 
While such SM processes introduce some degree of CP 
violation (CPV), the effects in the quark sector are far 
too weak to explain the observed matter dominance of 
the universe [7]. Consequently, it is important to search 
for further non-SM sources of CPV. 

Studies of neutral B meson oscillations, whereby a 
neutral meson changes into its own antiparticle via a 
box-diagram- mediated weak interaction [6], can provide 
a sensitive probe for such CPV processes. The semilep- 
tonic mixing asymmetry, defined as: 

q _ r(B° q -> B° q -> e+x) - v(B a q -» B° q -> e-x) 
° sl r(5° -> 5° -> i+x) + r(s° -> 50 -> e-x) ' ( ' 

allows the effects of any CP-violating processes to be di- 
rectly observed in terms of the resulting asymmetry of 
the decay products. Here i denotes a charged lepton of 
any flavor, and q represents the flavor of the non-6 valence 
quark of the meson. 

In the standard model, the semileptonic mixing asym- 
metry is related to the properties of the corresponding 
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B meson system, namely the mass difference AM, 
M(B° qH ) - M(B° qL ), the decay-width difference Ar g 
T{B° qL ) - r(B qH ), and the CP-violating phase <f> q , by: 



\M? 2 \ 



sin (f> q 



AM, 



-tan^, 



(2) 



agreement of the CP-violating phase <f> s with SM predic- 
tions. Given the current body of experimental evidence, 
improved measurements of both and afj are required 
in order to constrain the possible sources of new physics 
in B meson mixing and decay [23]. 

This article describes the measurement of the semilcp- 
tonic mixing asymmetry for B® mesons, 



Here the states B qH and B qL are the heavy and light 
mass eigenstates of the B meson system, which differ 
from the flavor eigenstates. Mf 2 and Y\ 2 are respectively 
the off-diagonal elements of the mass and decay matri- 
ces [6]. 

The standard model predictions [8] for both a s A and 
are very small: 



4 = 

T(B° 



B° 



£+DW-X)-T(B° 



B° 



(8) 



T(B° -^B a ^ (+D(*)-X) + T(B° 



B° 



i-D(*)+xy 



<4 = (-0.041 ±0.006)%, 
a s sl = (0.0019 ±0.0003)%. 



(3) 
(4) 



These predictions are effectively negligible compared to 
the current experimental precision. Hence, the measure- 
ment of any significant deviation from zero is an unam- 
biguous signal of new physics, which could lead to order- 
of-magnitude enhancements of \a^\ [9]. 

The B° semileptonic mixing asymmetry, a^, has been 
extensively studied by the B factories operating at 
the T(4S') resonance, including measurements by the 
CLEO [10, 11], BaBar [12, 13], and Belle [14] collabora- 
tions. The current world average of these measurements 
is [6]: 



without the use of initial-state flavor tagging. The flavor 
of the B° meson at the time of decay is determined by 
the charge of the muon in the semileptonic decay. Two 
separate decay channels are used: 

1. B° -> n+vD-X, 
with D~ -> K + ir~TT~ 

(plus charge conjugate process); 

2. B° — > n+isD*-X, 

with D*- -> D tt-,D q -> K+tt- 
(plus charge conjugate process); 



<4 = (-0.05 ±0.56)%. 



(5) 



Additional inclusive measurements from LEP [15-17] and 
DO [18] are subject to contamination from B° s mesons, 
and the extraction of relies upon assumptions about 



the contribution from afj. 



The two channels are treated separately, with each being 
used to extract a^, before the final measurements are 
combined. For clarity, the two channels are respectively 
denoted by fiD and \xD* throughout this paper, with the 
appropriate combinations of charges implied. Charges 
are only explicitly shown when required to describe the 
asymmetry measurement, or to avoid possible ambiguity. 



II. ANALYSIS OVERVIEW 

Experimentally, the semileptonic mixing asymmetry is 
expressed as: 



A- A 



BG 



pose 
r B° 



The recent evidence for a non-zero dimuon charge 
asymmetry by the DO experiment is sensitive to the lin- 
ear combination of B° and B® mixing asymmetries, with 
approximately equal contributions from each source [19]. 
The measurement constrains a band in the (<4,aJ) Here, A is the measured raw asymmetry, defined by: 
plane, which is inconsistent with the SM prediction at 
the 3.9 standard deviations level. By dividing the sam- 
ple into two components with different relative contribu- 
tions from B° and B®, the semileptonic asymmetries are 
measured to be: 



A = 



£>(*>- 



N, 



N, 



diff 



(9) 



(10) 



4(M = (-0.12 ±0.52)%, 
<i(w) = (-1.81 ± 1.06)%, 



(6) 
(7) 



where the measurements have a correlation coefficient 
of —0.799. The above extraction assumes that any new 
source of CPV entering the dimuon asymmetry does so 
through B mixing. Alternative hypotheses, for example 
new sources of dimuons from non-SM processes, cannot 
be excluded. 

Recent searches for CPV in B® — > J/ip(j) decays from 
the DO [20], CDF [21], and LHCb [22] collaborations find 



where N^± D (*) T is the number of reconstructed ^i ± I?(*) =F 
signal candidates. The sum is extracted by fitting the 
total mass distribution, and the difference by fitting the 
difference of two charge-specific mass distributions. The 
term ^4bg accounts for inherent detector-related back- 
ground asymmetries, for example due to the different 
reconstruction efficiencies for positively and negatively 
charged kaons. The denominator F^ S C is defined as the 
fraction of all fiD^ signal events that arise from decays 
of B° mesons after they have oscillated. It is required 
to account for mesons arising from direct B° de- 

cays, decays of B^- and B® mesons, or direct hadroniza- 
tion from cc quarks. All background asymmetries are 
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extracted using data-driven methods, while Monte Carlo 
(MC) simulation is used to determine the fraction of B° 
mesons that have undergone mixing prior to decay. 

This measurement assumes that the initial production 
of B°-B° is symmetric, and that there is no asymmetry 
in the decays of unmixed B° or B mesons (that would 
imply CPT violation), and no direct CP asymmetry in 
the semileptonic decays to charm states, or the decay 
of these charm states to the indicated products. With 
these assumptions, any observed semileptonic asymmetry 
would have to arise due to the mixing process. 

The B° meson has a mixing frequency AM,; = 0.507± 
0.004 ps _1 , of comparable scale to the lifetime t(B°) = 
1.518 ± 0.007 ps [6]. Hence the fraction of oscillated B° 
mesons is a strong function of the measured lifetime. The 
proper decay length ct for a particle is given by: 

L T cM T cM 
ct = -a-=L = L xy , 11 

h p pt 

where 7 and /? are the usual relativistic kinematic quan- 
tities; p, M and L are, respectively, the particle momen- 
tum, mass and decay length in the detector reference 
frame. The best precision is obtained by using the trans- 
verse quantities L xy and pt, due to finer instrumentation 
for tracking in this plane. The transverse decay length 
L xy is the projection of the vector pointing from the pro- 
duction to the decay vertex of the B meson onto the B 
meson transverse momentum direction. It can be nega- 
tive due to the limited spatial resolution of the detector. 

For semileptonic decays, the missing energy due to the 
undetected neutrino results in the measured transverse 
momentum being underestimated with respect to the ac- 
tual value. Hence the measured variable is actually the 
visible proper decay length (VPDL): 



VPDL(B) = L xy (B) 



cM(B) 
Pr{p.D)' 



(12) 



The dilution F^f is a very strong function of this vari- 
able, increasing monotonically with VPDL. To exploit 
this behavior, the measurements of all asymmetries are 
performed separately in bins of VPDL(_B°). These mea- 
surements are then combined for each channel to obtain 
the final measurement. The selected VPDL(_B°) bins are 
defined by the edges {-0.10, 0.00, 0.02, 0.05, 0.10, 0.20, 
0.60} cm. The /iD^ signal contributions outside of this 
range are found to be negligible. The first two bins in 
VPDL have negligible contributions from oscillated B° 
mesons, and are not included in the final measure- 
ment. They represent a control region in which the mea- 
sured raw asymmetry should be dominated by the back- 
ground contribution, i.e., A — Abg ~ 0. 

There can be significant (~1%) asymmetries due to 
detector effects. In particular, the material and detector 
elements that a particle traverses are different for pos- 
itively and negatively charged particles, as a result of 
the specific orientation of magnetic fields in the central 
tracking and muon detectors. In this analysis, such ef- 
fects are removed to first-order by reweighting all events, 



such that the total weight of events collected in each of 
the four (solenoid, toroid) magnet polarity configurations 
is the same (see Section III) . Remaining asymmetries are 
of order 0.1%, and are corrected using data-driven meth- 
ods. 

To avoid possible experimental bias, the central values 
of the raw asymmetries were hidden until all analysis 
methods were finalized. Initially this was achieved by 
randomly assigning all candidate charges; later, to allow 
the background asymmetries to be examined, the true 
charges were used, but unknown offsets were added to 
the raw charge asymmetries. 

The DO detector is briefly described in Section III, 
highlighting those features most relevant for this mea- 
surement. The event selection and raw asymmetry ex- 
traction are described in Sections IV- VI. The determi- 
nation of background asymmetries is described in Sec- 
tion VII, while Section VIII covers the extraction of the 
oscillated B° fraction. The results and conclusions are 
presented in Sections IX-XII. 



III. THE DO DETECTOR 

The DO detector has been described in detail else- 
where [24]. The most important detector components 
for this measurement are the central tracking system, the 
muon detectors, and the magnets. 

The central tracking system comprises a silicon mi- 
crostrip tracker (SMT) and a central fiber tracker (CFT), 
both located within a 2T superconducting solenoidal 
magnet. The SMT has «800,000 individual strips, with 
typical pitch of 50 — 80 /im, and a design optimized for 
tracking and vertexing capability at pseudorapidities of 
\r]\ < 2.5, where 77 = — ln[tan(0/2)] and 8 is the polar 
angle with respect to the beam axis. The system has a 
six-barrel longitudinal structure, each with a set of four 
layers arranged axially around the beam pipe, and in- 
terspersed with 16 radial disks. In the spring of 2006, 
a "Layer 0" barrel detector with 12288 additional strips 
was installed [25], and two radial disks were removed. 
This upgrade defines the chronological boundary between 
the two running periods, denoted Run Ha and Run lib. 
The sensors of Layer are located at a radius of 17 mm 
from the colliding beams. The CFT has eight thin coax- 
ial barrels, each supporting two doublets of overlapping 
scintillating fibers of 0.835 mm diameter, one doublet be- 
ing parallel to the collision axis, and the other alternating 
by ±3° relative to the axis. Light signals are transferred 
via clear fibers to solid-state photon counters that have 
«80% quantum efficiency. 

A muon system resides beyond the calorimeter, and 
consists of a layer of tracking detectors and scintillation 
trigger counters before a 1.8 T toroidal magnet, followed 
by two similar layers after the toroid. Tracking at \rj\ < 1 
relies on 10 cm wide drift tubes, while 1 cm mini-drift 
tubes are used at 1 < \rj\ < 2. 

The polarities of both the solenoidal and toroidal mag- 
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nets were regularly reversed during data acquisition, ap- 
proximately every two weeks, resulting in almost equal 
beam exposure in each of the four polarity configura- 
tions. This feature of the DO detector is crucial in reduc- 
ing detector-related asymmetries, for example due to the 
different trajectories of positive and negative muons as 
they traverse the magnetic fields in the detector. 



IV. EVENT SELECTION 

This analysis uses data collected by the DO detec- 
tor from 2002-2011, corresponding to approximately 
10.4 fb _1 of integrated luminosity, and representing the 
full Tevatron Run II sample of pp collisions at center-of- 
mass energy yfs = 1.96 TcV. Signal candidates are col- 
lected using single and dimuon triggers, and events that 
exclusively satisfy an impact-parameter based trigger are 
removed. 

For both channels, events are considered for selection if 
they contain a muon candidate with reconstructed track 
segments both inside and outside the toroid magnet. The 
muon candidate must be matched to a track in the central 
tracking system, with at least three hits in both the SMT 
and CFT. In addition, it must have transverse momen- 
tum px > 2 GeV/c, and total momentum p > 3 GeV/c. 

For events fulfilling these requirements, Z)(*)t candi- 
dates are constructed by combining three other tracks 
associated with the same initial pp interaction. Each 
track must satisfy pr > 0.7 GeV/c, and have at least two 
hits in both the SMT and CFT. The tracks must have 
a summed charge of magnitude \q\ = 1, with opposite 
sign to the muon charge. Each of the tracks comprising 
the like-charge pair is assigned the charged pion mass [6] . 
The third track, which has the same charge as the muon, 
is assigned the charged kaon mass [6] . 



A. fj,D Channel 

For the D~ — > K + tt~tt~ decay (and charge conjugate 
process) , the three hadron tracks must be consistent with 
originating at a single common vertex, with a vertex fit 
to the three tracks satisfying x 2 (vertex) < 16. These 
tracks are combined to construct & D~ candidate. The 
resulting D~ trajectory must be consistent with form- 
ing a common vertex with the muon to reconstruct a B° 
candidate. The cosine of the angle 9® between the mo- 
mentum and trajectory vectors of the D~ meson in the 
transverse plane must satisfy cos(#|?) > 0.0; i.e., the two 
vectors must point to the same hemisphere. The invari- 
ant masses must satisfy 1.6 < M(D~) < 2.1 GeV/c 2 and 
2.0 < M(B°) < 5.5 GeV/c 2 . 

At this preselection stage, a total of ~830 million can- 
didates remain, dominated by random three-track com- 
binations incorrectly associated with a real muon. A fit 
to the M(Knn) distribution yields 1 629 000 ± 29 000 \iD 



combinations. To increase the signal fraction of the sam- 
ple, a log likelihood ratio (LLR) method is utilized [26] 
to construct a single discriminating parameter from the 
combination of thirteen individual variables: the D~ 
transverse decay length L xy (D~), and its significance 
L xy (D~)/a[L xy (D~)]; the track isolation I of the kaon, 
the leading pion, and the trailing pion; the transverse mo- 
mentum of the kaon, the leading pion, and the trailing 
pion; the invariant mass of the reconstructed B candi- 
date, M(/j,D); the % 2 of the vertex fit for both the Ktttt 
and /j,D vertices; and the two-dimensional angular sepa- 
ration AR of the kaon and trailing pion, and of the two 
pions. The two-dimensional angular separation of two 
tracks is defined as AR — y/ 'A(f> 2 + Arf, where 77 is the 
pscudorapidity and <p is the azimuthal angle of each track. 
The track isolation / is the momentum of a particle di- 
vided by the sum of momenta of all tracks contained in a 
cone of size AR = 0.5 around the particle. Tracks corre- 
sponding to the other three final state particles for this 
candidate are excluded from the sum. 

The signal distributions required to construct the LLR 
discriminant are obtained from MC simulated events, in 
which the signal channel is required at generation, and 
the reconstructed tracks are required to match the correct 
particles at the generator level. For all MC studies de- 
scribed in this article, events are generated using PYTHIA 
version 6.409 [27], interfaced with evtgen [28] to model 
the decays of particles containing b and c quarks. Gen- 
erated events are processed by a GEANT based detector 
simulation, and overlaid with data from randomly col- 
lected bunch crossings to simulate pile-up from multiple 
interactions. The MC samples are then reconstructed us- 
ing the same software as used for data. The correspond- 
ing background distributions are obtained from sideband 
events in real data, defined by [1.660 < M(Ktttt) < 
1.760,1.964 < M(Kmr) < 2.064] GeV/c 2 , with each 
sideband scaled to give equal weight to the final distri- 
butions. 

Candidates enter the final data sample if the LLR dis- 
criminant exceeds a value £ m in, chosen to maximise the 
signal significance in data, Ns/\/Ns + Nb, where Ns 
and Nb are the number of /iD signal and background 
events, respectively. This figure of merit is found to cor- 
respond to the minimum uncertainty on the measured 
raw asymmetry. The optimal requirement is determined 
separately in each VPDL bin, and the value of L m - m de- 
creases for longer lifetimes, where the background from 
random track combinations is significantly reduced. 

After applying all selection requirements, the total fiD 
signal yield is ^740 000, with an overall efficiency of ap- 
proximately 44% with respect to the preselection sample. 
The signal efficiency in VPDL bins 3-6, used to extract 
afj, ranges from 53% to 72%. The M(Ktttt) invariant 
mass distribution over the full VPDL range is shown in 
Fig. 1. 
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FIG. 1: Distribution of the invariant mass M(K-nn) after all 
selections have been applied, for the [iD channel. The events 
have been weighted using the method described in Section V. 
The histogram shows the fit model used to extract the yield, 
with the background component drawn separately as a dashed 
line (see Section VI for fit models). 



B. /xD* Channel 

For the D*~ -> D°tt~ , D° — > K + ir~ decay (and charge 
conjugate process), D° candidates are reconstructed by 
combining a pair of oppositely charged tracks passing 
the criteria described above. The two tracks must form 
a common secondary vertex, and are used to reconstruct 
the parent D° candidate, which must satisfy pt(£>°) > 
0.7 GeV/c and \rj(D°)\ < 2.0, The invariant mass must 
lie in the range 1.7 < M(Ktt) < 2.0 GeV/c 2 . 

Next, an additional track is combined with the D° can- 
didate, which must have the opposite charge to the muon, 
and be consistent with forming a common vertex with the 
D°. This track is allocated the mass of the charged pion, 
and is here denoted 7Td*. The difference in the invari- 
ant masses of the D*~ and D° candidates must satisfy 
0.120 < [AM ee M(Kmr D ,)- M(Kir)\ < 0.200 GeV/c 2 . 
In addition, the displacement of the D° — > K + ir~ decay 
vertex with respect to the D*~ — > D°ir~ decay vertex 
must correspond to a significance of at least 3er, i.e., 



vW^t) 2 + (e L /a L ) 2 > 3, 



(13) 



where £t(L) an d ct(l) represent the distance and cor- 
responding uncertainty of the transverse (longitudinal) 
displacement between the two vertices. 

The D*~ candidate is then combined with the muon, 
to form a B° candidate. The muon, D , and 7td» tra- 
jectories must be consistent with arising from a common 
vertex, and the invariant mass of the B° must satisfy 
2.0 < MQiD*) < 5.5 GeV/c 2 . 

The final event selection requirement utilises a boosted 
decision tree (BDT) to further suppress backgrounds [29] . 
A total of 22 variables are selected as inputs: 



/(£>*), 7(5°); 

AR(K, n D * 



transverse 
Pt(D°); 



momentum px(K), Pt(^)i Pt(^d*), 



• angular separation AR(K, it) 
AR{ir,n D ,), AEpV); 

• transverse decay length L xy (D°), error a[L xy (D )], 
and significance L xy (D°)/a[L xy (D )}; 

• cosine of the angle, in the transverse plane, between 
the D° momentum vector and the position vector of 
the D decay vertex with respect to (a) the primary 
pp interaction vertex, and (b) the B° decay vertex; 

• cosine of the angle, in the transverse plane, between 
the D* momentum vector and the position vector 
of the D* decay vertex with respect to the primary 
vertex: 

• decay vertex fit quality x 2 (B°); and 

• invariant mass M(Kir) and M(/j,D*). 

The signal distributions are taken from MC simula- 
tion, in which the signal channel is generated exclusively 
by forcing the required decays in evtgen, and the recon- 
structed tracks are required to match the correct particles 
at the generator level. The background distributions are 
taken from real data, in which the kaon and two pions 
all have the same charge, and the muon has the opposite 
charge. The choice of BDT cut used to define the final 
data sample is made separately for Run Ha and Run lib 
samples, and for each VPDL(i?°) region, to optimise the 
signal significance in each case. 

After application of all selection criteria, the sample 
contains ^545 000 fiD* signal candidates. The AM dis- 
tribution for the full VPDL range is shown in Fig. 2. 



V. EVENT WEIGHTS 

In any given configuration of the solenoidal and 
toroidal magnet polarities, there can be detector-related 
asymmetries. These originate from differing detection 
efficiencies for positively and negatively charged parti- 
cles, in turn caused by their different trajectories as 
they bend through the magnetic fields in the detector. 
The regular reversal of both magnet polarities suppresses 
such effects. To ensure maximal cancellation of these 
instrumental asymmetries, an additional event-by-event 
weighting is applied such that the sums of weights in each 
(solenoid, toroid) configuration are the same, for a given 
sample. 

The weights are determined after applying the final 
event selections, by counting the total number of events 
in each of the four solenoid and toroid magnet polarity 
configurations. The weight for an event collected in a po- 
larity configuration i = {1, 2, 3, 4} is defined as N min /Ni 
where Ni is the number of events in this polarity configu- 
ration, and -/V m i n is the smallest of the four yields -^2,3.4- 
This procedure is performed separately for each channel, 
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FIG. 2: Distribution of the invariant mass difference AM = 
[A4(KmT£,*)~M(Kn)], for the \xD* channel. The events have 
been weighted using the method described in Section V. The 
solid line shows the fit model used to extract the yield. The 
triangular data points show the corresponding distribution 
for fiD* candidates in which the three hadrons have the same 
charge, scaled to give the same yield as the signal sample, 
in the sideband region 0.155 < AM < 0.170 GeV/c 2 (see 
Section VI for fit models). 



and for each VPDL(£?°) bin. Event weights are typically 
in the range 0.90-1.00, with very little variation between 
VPDL bins. For the unbinned samples, the total signal 
yields after event weighting are N(/j,D) = 721 519 ±3537, 
and N(^D*) = 519 066±3446, as shown in Figs. 1 and 2. 



VI. EXTRACTING THE RAW ASYMMETRY 

The raw asymmetry is extracted by fitting the invari- 
ant mass distributions M(Kttk) (or AM) for the £)W 
candidates. The sum distribution H sum is constructed by 
weighting all /j,D^*' candidates according to the magnet 
polarity weight. A difference distribution -f/diff is con- 
structed by taking the difference betwen the fi + D^*>~ 
and \jT £)(*)+ distributions. The sum and difference dis- 
tributions are modeled by, respectively, the functions: 



sum 



r sura — r sum T" Jv sum 



N sum ■ F si ^ (14) 
F diff = Ffog + A ■ A/sum • F sis , (15) 



where N sum is the total fj,D^ yield, and A is the corre- 
sponding raw charge asymmetry defined in Eq. (10). Dif- 
ferent models are used to parametrize the backgrounds 
for the sum (F s ^) and difference (-F^ff) histograms, 
while a single model F slg is used for the signal in both 
cases. The yields, asymmetries, and signal and back- 
ground parameters in these models are extracted by a 
simultaneous binned fit to the two distributions, to min- 
imise the total x 2 with respect to the fitting functions. 



A. fj,D Channel 

The physical width of the D~ meson is negligible com- 
pared to the detector resolution; therefore, the signal 
parametrization is chosen based on studies of simulated 
data to determine the mass resolution for this channel. 
The signal is modeled by the sum of two Gaussian func- 
tions constrained to have the same mean, but with dif- 
ferent widths and relative normalizations: 



F siB (D) 



1 



2tt 



1 



fa ' 



. ~(x-M D y/(2°i) 



(16) 



+ (1 - f G1 ) ■ — ■ e -{x-Mo) 2 /(2ol) 
CT2 



where x is the reconstructed invariant mass of the D~ 
candidate, Mn is the mean of the Gaussian peak, o\.2 
are the widths of the first and second Gaussians, and fci 
is the fraction of the signal in the first Gaussian peak. 

The background in the sum distribution exhibits 
slightly different behavior for each VPDL bin, hence a 
flexible parametrization is selected to provide good agree- 
ment in all bins, comprising the sum of three possible 
components: a low-order polynomial function, a falling 
exponential function, and a hyperbolic tangent function: 



(17) 



(N t , 



l)' \ /tanh ' -^tanh 

"H (1 ,/tanh) ' fpoly ' -^poly 

+ (1 - /tanh) ' (1 — /poly) ' ^exp}- 



Here /tanh arid / po iy are free parameters between zero 
and one, controling the relative contributions of the three 
components. 

The polynomial function includes a constant, linear, 
and cubic term: 

-Fpoly = Ci +pi ■ (X - Z m id) +P3 ' (x ~ X mid ) 3 , (18) 

where pi and P3 are free parameters of the fit, and a; m id is 
the mid-point of the fitting range. The exponential term 
is: 



F c: 



= G 2 • e m,n ', 



(19) 



where r is a free parameter and x m j n is the lower-limit of 
the fitting range. 

Finally, a hyperbolic tangent function described by 

^tanh = C 3 ■ [1 - tanh(fc • (x - M D ))\ (20) 

is used to model the effects of partially reconstructed de- 
cays, and reflections from decays into other three-track 
combinations. Monte Carlo simulations confirm that this 
parametrization is a good model for this source of back- 
ground, which includes contributions from D~ decays to 



-7T°, and K+K- 



four charged hadrons; and decays of D* 



D° decays to 
->■ D°tt-, with 



D" 



where the ir is not reconstructed. The 



9 



steepness of the threshold, denoted by k in Eq. (20), is 
controlled by the detector mass resolution. As such it is 
fixed according to the widths of the two Gaussian peaks: 

' (21) 



k = 



2 

mean 



where a mean = f G i ■ <J\ + (1 - /gi) ■ o" 2 is the weighted 
mean of the two widths. 

All three individual components are normalized to have 
unit area in the fitting range, by suitable choice of the 
constants C\, C2, and C3. The overall normalization 
is set by subtracting the fitted number of signal events 
(JV sum ) from the total event count in the sample (iVtot)- 
Using the total event count as a constraint in this way 
improves the fit precision. The free parameters are the 
two fractions, /tanh and / po iy, the two polynomial coef- 
ficients, £11,3, and the argument of the exponential func- 
tion, r. This empirical choice provides good agreement 
with the data, with relatively few free parameters, over 
a range of different background shapes. To improve fit 
precision and stability, each term in F s ^ is only used 
if it improves the fit probability. As a result, for VPDL 
bins 1-3 the exponential component is removed; for bin 
6, the cubic term is removed. 

For the difference fit, the overall normalization N^g of 
the background is fixed according to the observed number 
of D + and D~ events, and the signal contribution: 



A-N s . 



(22) 



where N^ t is the sum of all event weights for D ± candi- 
dates. The background model comprises the same three 
components as used in the fit to the sum. The thresh- 
old component is modeled by the same shape as in the 
sum fit, with the yield scaled by a free parameter (atanh) 
accounting for the possible charge asymmetry from this 
contribution. The combinatorial background is modeled 
by the sum of exponential and polynomial terms, with 
the shape parameters common to the sum fit, and the 
yield constrained by Eq. (22) after the contribution of 
the threshold component has been accounted for: 



(D) = 

N BG 



(23) 



diff ' { a tanh ' /tanh ' -Ftanh 

+ (I" Otanh ' /tanh) ' /poly ' -fpoly 



+ (1 a tanh ' /tanh) ' (1 /poly) ' -Fcxp}- 

This function has only one additional free parameter, 
with respect to the fit over the sum of all candidates, 
namely the asymmetry on the hyperbolic tangent, a tan h- 
The corresponding asymmetry term for the polynomial 
component is eliminated by applying the constraint from 
Eq. (22). 

In total, there are twelve free parameters in the mass 
fit for this channel, six describing the signal, and six de- 
scribing the background. The default fit is performed 
over the range [1.65 < M(Ktttt) < 2.05] GeV/c 2 , using 
100 bins of width 4 MeV/c 2 , with variations on both the 
fitting range and the bin width considered as sources of 
systematic uncertainty. 



B. (j,D* Channel 

For this channel the invariant mass difference distribu- 
tion AM = M(Kmr D , ) - M(Kn) is fitted to extract the 
raw asymmetry. The proximity to the pion production 
threshold at approximately 140 MeV/c 2 leads to phase- 
space effects that tend to distort the signal and back- 
ground distributions. To account for these effects, and 
based on studies of MC simulation data, the signal is 
modeled by a skewed triple-Gaussian function, i.e., three 
Gaussian peaks constrained to have the same mean, but 
with different widths and relative contributions, and each 
multiplied by a threshold shape: 



F sig (D*) = f G1 

+ (1 - /Gi) • f G 2 

+(1 - f G1 ) • (1 - / G2 ) 



G B k(x,ai,M,s) 
G sk {x, (T 2 , M, s) 
G sk (x,<r 37 M,s) 



(24) 



where G sk is the skewed Gaussian function: 



1 (t — M\ (s-M) 2 

G sk (x,a i ,M,s) = ^^-E[-s [ ■ ^ > }-e 



2-KGi 



V2<7 t 



(25) 



Here x is the reconstructed value of AM for this D*~ 
candidate, M is the mean of the Gaussian peak, <Ji, 2 ,3 
are the widths of the three Gaussians, and /gi,2 describe 
their relative fractional contributions. The function E(s- 
y) is a threshold shape modeled by the complementary 
error function, taking the skew s as an input parameter, 
and defined as: 



E(s ■ y) 



2 r 



dt. 



(26) 



The background in the sum fit is modeled by the prod- 
uct of a linear function and a power law function with a 
threshold at the charged pion mass M w [6] , and three free 
parameters a, b, and d: 



F^(D*) = d-{x-M„) a -(l + bx), 



(27) 



The background in the difference distribution is modeled 
by the same shape as used for the sum, but with a differ- 
ent overall scale, quantified by a background asymmetry 
parameter obg : 

F|f (£>•) = d ■ a BG ■ (x - ■ (1 + bx). (28) 



In total there are thirteen free parameters in the 
fit, nine for the signal, and four describing the back- 
ground. The default fit is performed over the range 
[0.139 < AM < 0.170] GeV/c 2 , using 62 bins of width 
0.5 MeV/c 2 , with variations on both the fitting range 
and the bin width considered as sources of systematic 
uncertainty. 
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C. Results 

The values of all physics parameters returned by the 
fits, for both channels, and for each of the six VPDL(£>°) 
bins are collected in Tables I— II, with examples of the fit 
projections shown in Fig. 3. Significant positive asym- 
metries are observed for all VPDL bins, including those 
in the control region VPDL(£?°) < 0.02 cm. This is ex- 
pected as a consequence of the positive kaon reconstruc- 
tion asymmetry, which is described and corrected for in 
Section VII. The two channels have similar statistical 
precision on the raw asymmetries, except for the first 
VPDL bin, where the sensitivity of the /j,D channel is 
significantly reduced by the increased background from 
random three track combinations close to the primary pp 
interaction. The \iD* channel is less susceptible to such 
effects, due to the intermediate resonance in the decay. 

The raw asymmetry measurement method is validated 
by the use of ensemble tests, in which the fits to data 
are repeated several thousand times with the \x^D^*^ 
charges randomized independently for each fit. Differ- 
ent input asymmetries are simulated, ranging from —5% 
to +5%, and the distribution of asymmetries extracted 
from the fits are examined. For all cases, the distribu- 
tions are well-modeled by Gaussian peaks, with a central 
value consistent with the input asymmetry, and a width 
consistent with the corresponding uncertainty reported 
in Tables I— II. A similar approach is used to confirm 
that the optimal precision is obtained by maximizing the 
signal significance. 



• for the nD* channel, the mass window used to se- 
lect D° — > Kn candidates is varied, from the nom- 
inal requirement of 1.7 < M(Ktt) < 2.0 GeV/c 2 , 
to twelve alternative ranges, giving different back- 
ground fractions and shapes in the final distribu- 
tion; 

• the function F slg , used to model the signal shape, is 
changed to an alternative choice; for the \iD case, a 
single Gaussian function is used; for the [iD* case, 
a skewed double-Gaussian is used; 

• the function F^g, used to model the background 
of the difference distribution, is changed to an al- 
ternative choice; for the \xD case, a linear function 
is used; for the fiD* case, a second or fourth order 
polynomial is used; 

• the event weights are allocated using an alternative 
method, based on the fitted number of fiD^ signal 
events in each polarity configuration, rather than 
the total number of candidates. 

To properly assess the combined effect of all these ad- 
justments to the fit, including correlations, all possible 
combinations of the above fit variations are tested, and 
the systematic uncertainty is allocated as the standard 
deviation of the full set of alternative measurements. Ta- 
ble III shows the final systematic uncertainties allocated 
for each VPDL bin, for both channels. The combined 
systematic uncertainty is significantly smaller than the 
statistical uncertainty in all cases. 



D. Systematic Uncertainties 

Several variations of the fits are performed to extract 
the raw asymmetry from the data, in order to examine 
the resulting spread of measured values, and assign an 
appropriate systematic uncertainty. These variations are 
considered for all cases where there is a reasonable alter- 
native to the choices made for the nominal fit, namely: 

• the lower and upper limits of the fitting region 
are varied to a number of different options within 
a 50 MeV/c 2 (10 McV/c 2 ) range of the nominal 
choice for the fj,D^ case; 

• the bin width is varied, with alternative widths of 
2-20 MeV/c 2 for the fiD case, and 0.5-2.0 MeV/c 2 
for the iiD* case; 

• for the fiD channel, the function Ff^ used to 
model the background of the sum distribution is 
changed to three alternative models in addition to 
the nominal choice; in one alternative, the polyno- 
mial is fixed to a linear function; in another, the 
polynomial is set to a quadratic function; in the 
third alternative, the mean of the hyperbolic tan- 
gent is allowed to be a free parameter, rather than 
being constrained to the mean D~ mass; 



VII. ACCOUNTING FOR DETECTOR 
ASYMMETRIES 

Both channels used in this measurement are recon- 
structed from the final state particles fj, ± K ± TT T TT T . In 
relating the measured raw asymmetry to the physical 
asymmetry under investigation, the effects of possible 
charge asymmetries in particle reconstruction must be 
considered. Neglecting asymmetries of second order or 
higher, the background asymmetry simplifies to: 

A BG = a» +a K - 2a v . (29) 

where the asymmetries a x are defined as the difference 
in reconstruction efficiency e for the positively and neg- 
atively charged particles: 



A. Kaon Asymmetry 

By far the largest background asymmetry to be taken 
into account is due to differences in the behavior of pos- 
itive and negative kaons as they traverse the detector. 
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TABLE I: Results ol the raw asymmetry fits for the /iD channel, in each of the six bins of VPDL(B°). The uncertainties are 
statistical, as returned by the fits. 
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1 A Q J_ n A 1 

1.4o zh U.41 


i on _i_ n qq 
l.zU zh U.oo 


l\/f(D\ (\Ra\T I n l \ 

ivl yu ) t^ivie V / C ) 


i oaa q _i_ n « 
louu.o zh U.O 


i occ o _i_ n o 
louO.o zt U.z 


i occ 7 _|_ n o 
louO. ( zh U.z 


loDO.y zh u.z 


i oat, A _|_ n 9 
loDO.4 zh U.z 


1 Qf: A q _i_ n q 
10D4.0 zh U.O 


(Tgi (MeV/c 2 ) 


21.6 ± 1.6 


19.0 ± 0.9 


18.6 ± 1.0 


18.3 ± 1.0 


18.0 ± 1.7 


15.0 ± 3.0 


(T G2 (MeV/c 2 ) 


39.6 ± 8.8 


32.5 ± 1.4 


30.1 ± 1.1 


29.6 ± 1.3 


28.6 ± 1.8 


27.0 ± 1.5 


fai 


0.67 ± 0.16 


0.417 ± 0.071 


0.361 ± 0.076 


0.363 ± 0.084 


0.33 ± 0.13 


0.16 ± 0.11 


OStanh (%) 


3.75 ± 2.86 


1.32 ± 1.04 


-0.02 ± 0.84 


-0.04 ± 0.79 


-0.24 ± 1.00 


-0.99 ± 3.66 


X7ndf 


172 / 190 


199 / 190 


218 / 190 


199 / 188 


213 / 188 


170 / 189 


X 2 (sum) /ndf 


81 / 92 


105 / 92 


114 / 92 


108 / 90 


123 / 90 


79 / 91 


X 2 (diff)/ndf 


91 / 98 


94 / 98 


104 / 98 


92 / 98 


90 / 98 


91 / 98 



TABLE II: Results of the raw asymmetry fits for the /iD* channel, in each of the six bins of VPDL(£?°). The uncertainties are 
statistical, as returned by the fits. 





Bin 1 


Bin 2 


Bin 3 


Bin 4 


Bin 5 


Bin 6 


VPDL(B U ) (cm) 


-0.10 - 0.00 


0.00 - 0.02 


0.02 - 0.05 


0.05 - 0.10 


0.10 - 0.20 


0.20 - 0.60 


N(pD") 


59 823 ± 2398 


151 585 ± 1677 


132 227 ± 1092 


104 463 ± 921 


58 409 ± 651 


12 029 ± 233 


A(%) 


1.82 ± 0.67 


1.10 ± 0.30 


0.94 ± 0.30 


1.38 ± 0.33 


2.11 ± 0.44 


0.55 ± 0.99 


M (MeV/c 2 ) 


145.08 ± 0.06 


145.07 ± 0.02 


145.03 ± 0.02 


144.99 ± 0.02 


144.97 ± 0.03 


145.00 ± 0.06 


(tgi (MeV/c 2 ) 


0.76 ± 0.11 


1.59 ± 0.12 


1.72 ± 0.09 


1.75 ± 0.12 


1.67 ± 0.15 


1.39 ± 0.23 


(JG2 (MeV/c 2 ) 


1.45 ± 0.30 


0.84 ± 0.03 


0.86 ± 0.02 


0.90 ± 0.03 


0.90 ± 0.04 


0.79 ± 0.11 


a G3 (MeV/c 2 ) 


3.78 ± 0.76 


3.45 ± 0.30 


3.89 ± 0.31 


3.87 ± 0.32 


3.74 ± 0.30 


2.92 ± 0.36 


fai 


0.32 ± 0.14 


0.39 ± 0.03 


0.41 ± 0.02 


0.38 ± 0.02 


0.38 ± 0.04 


0.44 ± 0.10 


SG2 


0.64 ± 0.05 


0.69 ± 0.05 


0.74 ± 0.04 


0.74 ± 0.04 


0.69 ± 0.06 


0.52 ± 0.17 


S 


0.368 ± 0.083 


0.470 ± 0.032 


0.510 ± 0.024 


0.541 ± 0.025 


0.574 ± 0.037 


0.508 ± 0.079 


a BG (%) 


1.25 ± 0.16 


1.17 ± 0.22 


-0.10 ± 0.39 


-0.19 ± 0.53 


0.18 ± 0.77 


-0.94 ± 1.41 


X 2 /ndf 


133 / 113 


138 / 113 


203 / 107 


153 / 107 


165 / 107 


150 / 107 


X 2 (sum) /ndf 


80 / 53 


84 / 53 


159 / 53 


94 / 53 


94 / 53 


77 / 53 


X 2 (diff)/ndf 


52 / 60 


54 / 60 


45 / 60 


59 / 60 


71 / 60 


73 / 60 



Negative kaons can interact with matter in the tracking 
system to produce hyperons, while there is no equive- 
lent interaction for positive kaons. As a result, the mean 
path length for positive kaons is longer, the reconstruc- 
tion efficiency is higher, and the kaon asymmetry a K is 
positive. 

The kaon asymmetry is measured using a dedicated 
sample of K*° — > K + tt~ decays, based on the tech- 
nique described in Ref. [30]. The K + tt~ and K~n + 
signal yields are extracted by fitting the charge-specific 
M(i4T ± 7r T ) distributions, and the asymmetry is deter- 
mined by dividing the difference by the sum. The track 
selection criteria are the same as those required for the 
signal channels, and all events must contain a 
muon passing the selections described in Section IV. 
Since the K*° channel includes a final state pion, of oppo- 
site charge to the kaon, the correction will also absorb any 



tracking asymmetry affecting pion reconstruction. One 
of the two a* terms in Eq. (29) is eliminated as a result. 

As expected, an overall positive kaon asymmetry is 
observed, of approximately 1% in this channel. A strong 
dependence on kaon momentum and absolute pseudora- 
pidity is found, and hence the final kaon asymmetry cor- 
rection to be applied in Eq. (29) is determined by the 
weighted average of a K {p(K), \rj(K)\] over the p(K) and 
\r)(K)\ distributions in the signal events: 

24 /N K \ 

i=l v ' 

where the sum is over eight bins in kaon momentum mul- 
tiplied by three bins in absolute pseudorapidity, and is 
the total number of signal candidates over all bins. The 
kaon asymmetry as a function of momentum is shown in 
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FIG. 3: Examples of the raw asymmetry fit for the two decays channels, for the fifth VPDL(_B°) bin corresponding to (0.10 < 
VPDL(B°) < 0.20) cm. The left plots show the sum distributions; the right plots show the difference distributions. In both 
cases, the solid line represents the total fit function, with the background part shown separately by the dashed line (see text). 



Fig. 4 for each of the three regions, and the val- 

ues and bin definitions are listed in Table IV. A relative 
systematic uncertainty of 5% is assigned to each bin to 
account for possible variations in the yield when different 
models are used to fit the signal and backgrounds in the 
A"* mass distribution. 

The kaon momentum distributions for each channel, 
within each |?7(AT)| region, and for each VPDL(i3 ) bin, 
are determined by fitting the appropriate invariant mass 
distribution in each of the eight kaon momentum bins, 
using the same parametrizations F sum as descibed in Sec- 
tion VI. Following studies over a range of fit variations, a 
relative systematic uncertainty of 3% (0.5%) is assigned 
on all [iD (fiD*) yields. 

The final corrections for each VPDL bin in both chan- 
nels are presented in Table V. The kaon corrections for 
the yuD* channel are slightly smaller than for the fiD 
channel, due to different kaon kinematics in the two de- 
cay topologies. 



B. Track Asymmetry 

Unlike kaons, positive and negative pions have almost 
identical interaction cross sections in matter. Any pos- 
sible asymmetry will be dominated by effects from track 
detection and reconstruction, which should be removed 
to first order by the magnet polarity weighting. 

The transverse momentum dependence of any residual 
tracking asymmetry is studied in Kg — > ir + ir~ decays. 
This channel can only be observed if a pair of oppositely 
charged pions is reconstructed, hence it is insensitive to 
the absolute asymmetry, and the overall scale is arbi- 
trarily fixed by setting the asymmetry in the lowest px 
bin to zero. The relative asymmetry as a function of 
Pt is determined by extracting the Kg yields in bins of 
[pt(7t + ),Pt(7I' _ )], and following the method described in 
Ref. [30], except that K s — > n + n~ decays are used in- 
stead oij/ip —> n + fJ.~ ■ As shown in Fig. 5, no evidence of 
any track px dependence is observed, over the range 0.5- 
7 GeV/c, within an uncertainty of ±0.05%. As a result, 
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TABLE III: Systematic uncertainties on the raw asymmetry measurement for both channels, extracted from an ensemble of 
fits with variations on each quantity tested. Also shown are the corresponding statistical uncertainties, for comparison. 



Source 


Bin 1 


Bin 2 


Bin 3 


Bin 4 


Bin 5 


Bin 6 


-0.10 - 0.00 cm 


0.00 - 0.02 cm 


0.02 - 0.05 cm 


0.05 - 0.10 cm 


0.10 - 0.20 cm 


0.20 - 0.60 cm 


fiD channel 


Bin width 


0.09% 


0.01% 


0.01% 


0.01% 


0.00% 


0.05% 


Fit limits 


0.17% 


0.06% 


0.08% 


0.05% 


0.03% 


0.12% 


Magnet weighting 


0.02% 


0.00% 


0.00% 


0.00% 


0.00% 


0.01% 


Signal model 


0.03% 


0.03% 


0.01% 


0.04% 


0.01% 


0.01% 


Background model (sum) 


0.03% 


0.00% 


0.01% 


0.01% 


0.01% 


0.00% 


Background model (diff) 


0.01% 


0.00% 


0.01% 


0.00% 


0.01% 


0.02% 


Combined systematic 


±0.19% 


±0.07% 


±0.08% 


±0.07% 


±0.05% 


±0.13% 


Statistical 


±1.28% 


±0.35% 


±0.32% 


±0.33% 


±0.41% 


±0.88% 








channel 








Bin width 


0.06% 


0.03% 


0.02% 


0.02% 


0.01% 


0.08% 


Fit limits 


0.05% 


0.01% 


0.01% 


0.01% 


0.06% 


0.06% 


Magnet weighting 


0.01% 


0.01% 


0.00% 


0.00% 


0.00% 


0.01% 


Signal model 


0.01% 


0.01% 


0.05% 


0.06% 


0.03% 


0.03% 


Background model 


0.13% 


0.01% 


0.00% 


0.01% 


0.07% 


0.01% 


M(D°) cut 


0.01% 


0.01% 


0.01% 


0.01% 


0.02% 


0.02% 


Combined systematic 


±0.13% 


±0.04% 


±0.05% 


±0.07% 


±0.08% 


±0.09% 


Statistical 


±0.67% 


±0.30% 


±0.30% 


±0.33% 


±0.44% 


±0.99% 
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FIG. 4: Kaon asymmetry as a function of kaon momentum, 
for three regions of absolute pseudorapidity, as extracted from 
the K*° -> K + ir~ channel. 



TABLE IV: Kaon charge asymmetries in bins of p(K), as 
extracted from the K*° —¥ K + tv~ channel, for each of the 
three regions in absolute pseudorapidity: central < 
0.7); mid-range (0.7 < \v{K)\ < 1.2); and forward (1.2 < 
\V(K)\ < 2.2). 



p(K) 
range 


a K (%) 


Central 


Mid- Range 


Forward 


0.7 - 1.7 


1.38 ± 0.11 


1.38 ± 0.18 


1.23 ± 0.49 


1.7 - 2.4 


1.09 ± 0.14 


1.03 ± 0.16 


1.47 ± 0.23 


2.4 - 3.2 


0.76 ± 0.15 


0.78 ± 0.16 


1.53 ± 0.20 


3.2 - 4.2 


0.65 ± 0.18 


0.80 ± 0.18 


1.52 ± 0.20 


4.2 - 5.5 


0.78 ± 0.21 


1.36 ± 0.20 


1.51 ± 0.22 


5.5 - 7.5 


0.50 ± 0.22 


1.45 ± 0.20 


1.72 ± 0.21 


7.5 - 11.5 


0.24 ± 0.28 


1.41 ± 0.22 


1.16 ± 0.20 


> 11.5 


0.64 ± 0.38 


1.33 ± 0.26 


0.89 ± 0.20 



any residual tracking asymmetry will cancel to first order 
in the reconstruction of the pion and oppositely-charged 
muon, which remain to be taken into account after apply- 
ing the kaon asymmetry correction. There are insufficient 
statistics in the Kg — » it + tt~ channel to extend to higher 
transverse momenta. However, this momentum region 
contains the majority of fiD'*) signal candidates. 

A second dedicated channel K ± * — > K^ir^ is used 
to measure the absolute residual track asymmetry. The 
Kgir^ yields for each pion charge are extracted by fit- 
ting the Af(i^g7r ± ) invariant mass distributions, and the 



asymmetry calculated from the sum and difference of 
these yields. No significant asymmetry is found in this 
study, which is consistent with the findings of previous 
studies [30]. As such, no correction is applied to the 
asymmetry to account for the effects of track reconstruc- 
tion, and a lr in Eq. (29) is assigned to be zero. We allo- 
cate a systematic uncertainty of ±0.05% to account for 
the limited precision of this asymmetry measurement. 
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6 8 10 

P T (7t±) (GeV/c) 

FIG. 5: Relative track reconstruction asymmetry as a func- 
tion of track transverse momentum, as measured from Kg — > 
■k + -k~ decays. The absolute scale is chosen to give zero asym- 
metry for the first bin, since this channel is only sensitive to 
the variations in asymmetry between bins. 




P T (H±) (GeV/c) 



FIG. 6: Muon asymmetry as a function of muon transverse 
momentum, as extracted from a dedicated sample of J/tp — > 
H + jjT decays. 



C. Muon Asymmetry 

The residual charge asymmetry for muon identifica- 
tion is measured using J ftp — > /i + fi~ decays, using the 
technique developed in Ref. [30]. A small but signifi- 
cant asymmetry is observed, with a sizeable dependence 
on the muon transverse momentum, as shown in Fig. 6. 
The corresponding correction a M to be applied to the raw 
asymmetry is extracted using the same method as for 
the kaon asymmetry, by performing a weighted average 
of the muon asymmetry over bins of px (m) , analogous to 
Eq. (31). 

The final muon asymmetry corrections for each VPDL 
bin and both channels are summarized in Table V. To 
account for possible systematic uncertainties, the entire 
procedure of measuring muon asymmetries and convo- 
luting with the transverse momentum distributions is re- 
peated with several variations to the method, and the 



corresponding changes in the final measured muon asym- 
metry in each bin are used to assign a systematic uncer- 
tainty. The variations include changing the mass bin- 
ning of the M(fj, + distributions, changing the fitting 
function used to extract the J/rp yields, changing the px 
binning scheme, including an absolute pseudorapidity de- 
pendence, and using an alternative method to determine 
the polarity-based event weights. 



VIII. SAMPLE COMPOSITION: DILUTION 
FROM SYMMETRIC PROCESSES 

Not all /j,D^ combinations originate from the decay 
of oscillated B° mesons. Alternative charge symmetric 
sources will contribute only to the denominator in the 
raw asymmetry extraction, and hence dilute any physical 
asymmetry afj. 

In general, mesons containing a charm quark can be 
produced in many different ways, which we divide into 
five categories for the purposes of this measurement: 

1. direct hadroniszation from an initial c(c) quark, 
here denoted as 'prompt'; 

2. as a product of B° meson decay; 

3. as a product of B ± meson decay; 

4. as a product of B® meson decay; and 

5. as a product of a b baryon decay. 

The contribution from b baryons is found to be negli- 
gible, using the technique described below, and will no 
longer be considered. This scheme includes possible in- 
termediate excited resonances of both B and D mesons, 
for example processes such as c — > D*°X —> D + X' or 
those with higher excitations. For both neutral B meson 
sources, there may be mixing via box diagrams prior to 
decay, so sources 2 and 4 in the above list can be subdi- 
vided into 'mixed' and 'direct' decays. 

The total fraction of signal events coming from B° me- 
son decays is determined using inclusive MC simulations 
in which the only requirement at the generator level is 
the presence of the appropriate decay channel, and 

the presence of a muon (of any charge). The samples in- 
clude non-primary gluon splitting into heavy flavor bb 
and cc pairs, in addition to pair production from flavor 
excitation and flavor creation mechanisms. The gener- 
ated events are passed through the full simulation chain, 
and then processed by the same reconstruction and se- 
lection algorithms as used to select events from real data 
for the two signal channels. 

At the reconstruction level, the final state tracks must 
correspond to the true kaons and pions from the D^*^ 
decay. The result for each channel is a sample of 
events, with an accompanying reconstructed muon, in 
which the decay chain can be investigated in detail to 
extract the parentage information. The candidates are 
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TABLE V: The background asymmetries from kaon and muon reconstruction, in bins of VPDL(B°), for both signal channels. 
Also shown are the combined background asymmetries a K + a M , and the raw asymmetries. In each case, the upper uncertainty 
is statistical, the lower systematic. 





Bin 1 
-0.10 - 0.00 cm 


Bin 2 
0.00 - 0.02 cm 


Bin 3 
0.02 - 0.05 cm 


Bin 4 
0.05 - 0.10 cm 


Bin 5 
0.10 - 0.20 cm 


Bin 6 
0.20 - 0.60 cm 


/j,D channel 


A(%) 


2.70 ± 1.28 
± 0.19 


1.02 ± 0.35 
± 0.07 


1.16 ± 0.32 
± 0.08 


1.50 ± 0.33 
± 0.07 


1.48 ± 0.41 
± 0.05 


1.20 ± 0.88 
± 0.13 


a K (%) 


1.128 ± 0.041 
± 0.014 


1.124 ± 0.040 
± 0.014 


1.141 ± 0.040 
± 0.014 


1.147 ± 0.040 
± 0.014 


1.157 ± 0.040 
± 0.015 


1.157 ± 0.040 
± 0.014 


a" (%) 


0.102 ± 0.025 
± 0.008 


0.105 ± 0.027 
± 0.009 


0.107 ± 0.029 
± 0.012 


0.107 ± 0.029 
± 0.013 


0.108 ± 0.028 
± 0.011 


0.108 ± 0.028 
± 0.009 


a K + a" (%) 


1.230 ± 0.048 
± 0.016 


1.229 ± 0.048 
± 0.017 


1.248 ± 0.049 
± 0.018 


1.254 ± 0.049 
± 0.019 


1.265 ± 0.049 
± 0.019 


1.265 ± 0.049 
± 0.017 


fiD* channel 


A(%) 


1.82 ± 0.67 
± 0.13 


1.10 ± 0.30 
± 0.04 


0.94 ± 0.30 
± 0.05 


1.38 ± 0.33 
± 0.07 


2.11 ± 0.44 
± 0.08 


0.55 ± 0.99 
± 0.09 


a K (%) 


1.089 ± 0.047 
± 0.013 


1.078 ± 0.052 
± 0.014 


1.078 ± 0.050 
± 0.014 


1.085 ± 0.050 
± 0.014 


1.086 ± 0.049 
± 0.014 


1.098 ± 0.050 
± 0.014 


a" (%) 


0.097 ± 0.027 
± 0.012 


0.098 ± 0.031 
± 0.022 


0.101 ± 0.033 
± 0.023 


0.101 ± 0.033 
± 0.022 


0.101 ± 0.033 
± 0.020 


0.101 ± 0.031 
± 0.016 


a K + a" (%) 


1.186 ± 0.054 
± 0.018 


1.176 ± 0.061 
± 0.026 


1.179 ± 0.060 
± 0.027 


1.186 ± 0.060 
± 0.026 


1.187 ± 0.059 
± 0.024 


1.199 ± 0.059 
± 0.021 



weighted according to their true decay time, to ensure 
that the B meson lifetimes match the current world av- 
erage measurements, with the uncertainties on these life- 
times taken into account when assigning systematic un- 
certainties. 

Table VI lists the resulting fractions of fiD^ candi- 
dates from each source in both signal channels, as a 
function of the reconstructed VPDL(£>°). In general, 
the B° fraction is approximately 80-90%, except in the 
first (negative VPDL) bin in which the prompt contribu- 
tion reduces this to around 60%. The B ± contribution 
is small but significant, building from approximately 5- 
18% (6-14%) in the fiD^ case as the VPDL increases. 
This graduation is due to the longer lifetime of the B ± 
meson relative to B°. 

For all bins, the B° s fraction is very small, at approx- 
imately 1-3%. We correct for the possible contribution 
from the scmilcptonic mixing asymmetry in B® mesons, 
a|j, on a VPDL bin- by-bin basis, by extending Eq. (9) to 
include the non-zero B® fraction: 



A- A 



BG 



rposc 



pose 
r B° 



(32) 



Here F° S C = F(B° S ->• fiD^X) ■ Xs is the fractional con- 
tribution of oscillated B® mesons in the sample, where 
Xs = 0.499292 ± 0.000016 is the integrated mixing prob- 
ability [6]. The parameter is assigned the world aver- 
age value, [—1.05 ± 0.64]% [6]. The uncertainty on this 
quantity is taken into account when determining the sys- 



tematic uncertainty on the final measurement. 

The fraction of B° mesons that oscillate into their 
antiparticle prior to decay is determined by applying a 
weight Wp lK to all B° -> [iD^ events based on the true 
decay time (tj) of the B° meson: 



W 



, =-[l-cos(AM d ■*,•)], 



(33) 



where AMj is the mass difference of the heavy and light 
eigenstates in the B° system, assigned to be the world- 
average value 0.507 ± 0.004 ps _1 [6], with the precision 
taken into account when assigning systematic uncertain- 
ties. The resulting fractions Fg s o, are defined as the sum 
of these mixing weights divided by the total number of 
events in the MC sample. Table VI lists the resulting 
fractions versus VPDL for both channels. 

Various sources of systematic uncertainty on F^q are 
considered. The prompt fraction is negligible in VPDL 
bins 3-6 used in the final measurement; therefore, no 
systematic uncertainties are allocated from this source. 

The simulation may not describe the data perfectly 
In particular, the simulation doesn't account for any ef- 
fects due to the muon triggers used to collect data. MC 
simulations show that the pre-trigger muon transverse 
momenta distributions from B° and decays are com- 
pletely consistent, as expected from the closeness of the 
meson masses. On the other hand, there are small differ- 
ences in the Pt{ij) distributions for B® decays. Reweight- 
ing events by a trigger acceptance correction leads to a 
small reduction in the B° fraction, of order 3%. Since 
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TABLE VI: Fraction of L> w± candidates arising from each source as determined from MC simulation, for both channels, and 
in bins of VPDL(B°). For each channel, the upper row lists the central values and statistical uncertainties (from limited sample 
size in simulation), while the lower row lists the systematic uncertainties, only determined for the final F b s q values. 





Bin 1 
-0.10 - 0.00 cm 


Bin 2 
0.00 - 0.02 cm 


Bin 3 
0.02 - 0.05 cm 


Bin 4 
0.05 - 0.10 cm 


Bin 5 
0.10 - 0.20 cm 


Bin 6 
0.20 - 0.60 cm 


jiD channel 


F{c -> nDX) 
F(B° -> fiDX) 
F{B ± -> nDX) 
F(B° -> nDX) 


0.361 ± 0.011 
0.019 ± 0.003 
0.052 ± 0.005 
0.569 ± 0.011 


0.069 ± 0.003 
0.019 ± 0.001 
0.075 ± 0.003 
0.837 ± 0.004 


0.003 ± 0.000 
0.029 ± 0.001 
0.101 ± 0.003 
0.868 ± 0.003 


0.000 ± 0.000 
0.027 ± 0.001 
0.118 ± 0.003 
0.854 ± 0.003 


0.000 ± 0.000 
0.030 ± 0.002 
0.141 ± 0.004 
0.829 ± 0.004 


0.000 ± 0.000 
0.032 ± 0.004 
0.186 ± 0.008 
0.781 ± 0.009 


T7IOSC 


0.018 ± 0.003 
± 0.001 


0.009 ± 0.001 
± 0.000 


0.057 ± 0.002 
± 0.001 


0.208 ± 0.003 
± 0.005 


0.520 ± 0.005 
± 0.011 


0.658 ± 0.010 
± 0.017 


fiD* channel 


F(c^ nD*X) 
F(B° -> iiD'X) 
F(B ± -> ^D*X) 
F(B° -> fj,D*X) 


0.373 ± 0.010 
0.009 ± 0.002 
0.058 ± 0.005 
0.560 ± 0.010 


0.082 ± 0.003 
0.011 ± 0.001 
0.073 ± 0.003 
0.835 ± 0.004 


0.005 ± 0.001 
0.014 ± 0.001 
0.080 ± 0.003 
0.901 ± 0.003 


0.000 ± 0.000 
0.013 ± 0.001 
0.083 ± 0.003 
0.903 ± 0.004 


0.001 ± 0.000 
0.016 ± 0.002 
0.104 ± 0.005 
0.880 ± 0.005 


0.000 ± 0.000 
0.017 ± 0.005 
0.146 ± 0.013 
0.836 ± 0.013 




0.013 ± 0.002 
± 0.001 


0.010 ± 0.001 
± 0.000 


0.061 ± 0.003 
± 0.002 


0.231 ± 0.005 
± 0.003 


0.570 ± 0.008 
± 0.007 


0.713 ± 0.016 
± 0.009 



this source accounts for less than 3% of all candi- 
dates, the effect of this trigger correction on Fg„ is tiny, 
of order 0.001, and is neglected. 

The decay branching ratios of B° mesons into semilep- 
tonic final states containing a D T (D* T ) meson are 
known to around 10% (5%) precision [6]. As such, we 
vary the B° fractions up and down for the two channels 
by these fractions and assign a systematic uncertainty 
from this source equal to the total variation with respect 
to the default value. 

While the simulation assumes equal production frac- 
tions for the B° and mesons, based on isospin in- 
variance, experimental evidence from T(4S I ) decays can 
only confirm this assumption with a precision of around 
1% [6]. As such, we vary the B^ fraction up and down 
by this amount (and the converse variation for the B° 
fraction), and take the resulting change in Fgo c as a sys- 
tematic uncertainty. 

To account for the uncertainties on the world-average 
B meson lifetime values, we repeat the evaluation of F^q 
with the input lifetimes adjusted within their uncertain- 
ties, and assign a systematic uncertainty equal to the 
maximum deviation from the nominal Fg s „ value. Sim- 
ilarly, systematic uncertainties are allocated to account 
for the limited precision of AMd and aj. The break- 
down of systematic uncertainties is shown in Table VII. 
For the final measurement of a^, the uncertainties from 
the limited sample size in simulation are also categorized 
as systematic, not statistical, since they are not related 
to the size of the data sample. 



IX. RESULTS 

From the raw asymmetries and detector-related asym- 
metries listed in Table V, and the corresponding dilu- 
tion fractions Fg s o presented in Table VI, the final value 
of the semileptonic mixing asymmetry is determined 
for each VPDL bin i = {3, 4, 5, 6}, and for each channel 
3 = {fi>D,fj,D*}: 

_ A(ij) - a K (ij) - a»(ij) - F^jzj) ■ < 

(34) 

The first two VPDL bins, i — {1,2} are not included, as 
these represent the control region in which the expected 
signal contribution is negligible. 

To extract the corresponding uncertainties (both sta- 
tistical and systematic) care must be taken to properly 
account for all the correlations and constraints on the 
various inputs to the measurement. In particular: 

• the raw asymmetries for each bin and channel are 
independent; 

• the kaon asymmetry as a function of [p(K ),\-q(K)\] 
(Fig. 4) is 100% correlated between bins, and be- 
tween channels; 

• the muon asymmetry as a function of Pt(aO (Fig- 6) 
is 100% correlated between bins, and between chan- 
nels; 

• the asymmetry a^, used to derive the correction for 
a possible contribution from B® mixing, is 100% 
correlated between bins, and between channels; 
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TABLE VII: Systematic uncertainties from different sources on the dilution fraction Fgo, for both channels, and in bins of 
VPDL(B°). 





Bin 1 
-0.10 - 0.00 cm 


Bin 2 
0.00 - 0.02 cm 


Bin 3 
0.02 - 0.05 cm 


Bin 4 
0.05 - 0.10 cm 


Bin 5 
0.10 - 0.20 cm 


Bin 6 
0.20 - 0.60 cm 




Branching Ratios 
Production Fractions 
B meson lifetimes 
AM d 


± 0.001 
± 0.000 
± 0.000 
± 0.000 


± 0.000 
± 0.000 
± 0.000 
± 0.000 


± 0.001 
± 0.000 
± 0.000 
± 0.001 


± 0.004 
± 0.001 
± 0.001 
± 0.003 


± 0.009 
± 0.002 
± 0.003 
± 0.005 


± 0.015 
± 0.003 
± 0.007 
± 0.002 


Total 


± 0.001 


± 0.000 


± 0.001 


± 0.005 


± 0.011 


± 0.017 




Branching Ratios 
Production Fractions 
B meson lifetimes 
AM d 


± 0.001 
± 0.000 
± 0.000 
± 0.000 


± 0.000 
± 0.000 
± 0.000 
± 0.000 


± 0.001 
± 0.001 
± 0.001 
± 0.001 


± 0.001 
± 0.001 
± 0.001 
± 0.003 


± 0.004 
± 0.001 
± 0.003 
± 0.005 


± 0.006 
± 0.002 
± 0.005 
± 0.003 


Total 


± 0.001 


± 0.000 


± 0.002 


± 0.003 


± 0.007 


± 0.009 



• the oscillation fractions, F^„, are treated as inde- 
pendent. 

To ensure that all such correlations are taken into ac- 
count, the final statistical and systematic uncertainties 
on each measurement, and on the combination, are 
derived from 200 000 ensemble tests in which all input 
variables are randomly chosen according to a Gaussian 
probability density function, with an appropriate central 
value and width, and the distributions of the resulting 
dgj measurements are inspected and fitted. This process 
is performed twice, once with the inputs varied accord- 
ing to their statistical uncertainties, and once with the 
inputs varied according to their systematic uncertainties. 

Figure 7 and Table VIII show the individual results for 
the four signal VPDL bins in each channel, with statisti- 
cal and systematic uncertainties. 

Once the uncertainties on the individual measure- 
ments are established, the combination between VPDL 
bins, and then between channels, is performed. For each 
channel, the combined ajj value is obtained by a weighted 
average of the four individual measurements: 



„d U s _ Ef=3 a rffa>fa") 

a siu; - ^ 7~~, ■ 



(35) 



where the weights w{ij) are the inverse of the sum in 
quadrature of statistical and systematic uncertainties for 
that measurement: 



w(ij) 



1 



(36) 



The central values and uncertainties for the combinations 
are again determined by performing the full set of 200,000 
ensemble tests, with all inputs varied, and examining the 
effect on the final values of £jjj from each channel. This 



procedure yields the following results: 

a^(fiD) = [0.43 ± 0.63 (stat.) ± 0.16 (syst.)]%, (37) 
a^(fiD*) = [0.92 ±0.62 (stat.) ±0.16 (syst.)]%. (38) 

Finally, the combination is extended to give the full 
weighted average of the two channel-specific measure- 
ments, with full propagation of uncertainties, to yield 
the final measurement: 

4 = [0.68 ±0.45 (stat.) ±0.14 (syst.)]%. (39) 

The weights w(ij) used for this combination are pre- 
sented in Table VIII. 



X. CROSS-CHECKS 

To test the robustness of the measurement technique, 
the analysis is repeated with the event samples divided 
into pairs of orthogonal sub-sets, of approximately equal 
size. The raw asymmetries A, detector-related back- 
ground corrections a K and a M , and oscillation fractions 
.Fgo c are redetermined for each sub-sample, and the 
semilcptonic mixing asymmetry measured in each 
case. The sub-samples are defined by the following crite- 
ria: 

• r](fi) < and r](p) > 0; 

• \r]{K)\ < 0.7 and \r]{K)\ > 0.7; 

• p(K) < 3.2 GeV/c and p(K) > 3.2 GcV/c; 

• a chronological division corresponding to early and 
late data collection; 

• cr(VPDL) < 40 urn and cr(VPDL) > 40 /jm. 
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TABLE VIII: Individual measurements of (A — Abg) and af] in each of the six VPDL(_B°) bins, for both channels in this 
analysis. For each entry, the upper uncertainty is statistical, the lower systematic, with all correlations taken into account. 
Also shown are the weights iv(ij) used to combine the eight separate measurements, normalized to unity. 





Bin 1 


Bin 2 


Bin 3 


Bin 4 


Bin 5 


Bin 6 




-0.10 - 0.00 cm 


0.00 - 0.02 cm 


0.02 - 0.05 cm 


0.05 - 0.10 cm 


0.10 - 0.20 cm 


0.20 - 0.60 cm 


/j,D Channel 


A - Abg (%) 


1.48 ± 1.28 
± 0.20 


-0.20 ± 0.35 
± 0.09 


-0.07 ± 0.32 
± 0.10 


0.26 ± 0.33 
± 0.09 


0.23 ± 0.41 
± 0.07 


-0.05 ± 0.89 
± 0.14 


4 (%) 






-1.29 ± 5.68 


1.25 ± 1.61 


0.44 ± 0.79 


-0.07 ± 1.36 


Not used 


± 1.69 


± 0.43 


± 0.14 


± 0.21 


weight w(ij)/ J2ij 






0.006 


0.072 


0.309 


0.105 


fiD* Channel 


A - Abg (%) 


0.64 ± 0.67 
± 0.14 


-0.07 ± 0.31 
± 0.07 


-0.23 ± 0.31 
± 0.08 


0.20 ± 0.34 
± 0.09 


0.93 ± 0.44 
± 0.10 


-0.63 ± 0.99 
± 0.11 


< {%) 


Not used 


-3.79 ± 5.00 
± 1.27 


0.87 ± 1.45 
± 0.39 


1.63 ± 0.78 
± 0.17 


-0.89 ± 1.39 
± 0.15 


weight ■w(ij)/J2a w ( i j) 






0.007 


0.088 


0.311 


0.102 
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FIG. 7: Final measurements of the semileptonic asymmetry a^, in bins of VPDL(B°), for both channels. The cross-hatched 
bands show the mean values (and their total uncertainties) determined for each channel separately. 



The results are summarized in Table IX. In all cases, the 
measured values of are statistically consistent with 
each other, despite some samples having significantly dif- 
ferent background corrections. 

In addition, the measurement is repeated using only 
events that satisfy a single muon trigger. This corre- 
sponds to approximately 90% of the total sample. The 
resulting ai value for these events is consistent with the 
nominal value, taking into account the correlation be- 
tween the samples. 

The fraction of events from mixed B° decays, F^ s o, 
is a strong function of the visible proper decay length 
of the reconstructed B° candidate. Hence any non-zero 
value of afj will lead to a VPDL dependence on the back- 
ground subtracted asymmetry (A — Abg) [see Eq. (9)]. 



Figure 8 shows this dependence for both channels, with 
the Fg s o ■ aj~] distribution superimposed on the plot for 
comparison, using the final cq measurement from the two 
channel combination. The two distributions are statisti- 
cally consistent, indicating that the VPDL dependence of 
the observed background-subtracted asymmetry is con- 
sistent with the hypothesis that it originates from the 
mixing of B° mesons. The \ 2 quantifying this agreement 
between the (A — Abg) and Fg s a ■ distributions is 2.3 
(4.5) for the /j,D^ channel, compared to 2.7 (6.9) under 
the SM assumption for a^. For this test, the statistical 
and systematic uncertainties are combined in quadrature. 

The same data can be used to validate the en mea- 
surement using an alternative method, in which the dis- 
tribution of (A - Abg) versus VPDL(B°) is fitted to the 
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TABLE IX: Results of the analysis cross-checks, in which the 
data are divided into pairs of orthogonal and independent 
subsets, and the measurements of repeated for each sam- 
ple. The uncertainties shown here are the sum in quadrature 
of the statistical and systematic components. 



Sub-sample 
requirement 


4 (%) 


/iD channel 


/iD* channel 


Comb. 


Nominal Result 


0.43 ± 0.65 


0.92 ± 0.64 


0.68 ± 0.47 


< o 
v(p) > o 


0.38 ± 0.88 
0.53 ± 0.91 


0.60 ± 0.88 
1.21 ± 0.88 


0.49 ± 0.63 
0.88 ± 0.64 


\V(K)\ < 0.7 
\ V (K)\ > 0.7 


0.48 ± 0.95 
0.36 ± 0.85 


-0.39 ± 1.14 
1.17 ± 0.86 


0.04 ± 0.77 
0.76 ± 0.62 


p(K ) < 3.2 GeV/c 
p(K) > 3.2 GeV/c 


0.02 ± 0.87 
1.11 ± 0.92 


-0.30 ± 1.42 
1.00 ± 0.79 


-0.14 ± 0.85 
1.05 ± 0.62 


cr(VPDL) < 40 fim 
cr(VPDL) > 40 trni 


0.22 ± 0.84 
0.76 ± 0.95 


0.18 ± 0.95 
0.98 ± 1.01 


0.20 ± 0.65 
0.87 ± 0.70 


First half data 
Second half data 


0.82 ± 0.89 
0.19 ± 0.86 


1.39 ± 0.88 
-0.30 ± 1.02 


1.11 ± 0.67 
-0.06 ± 0.68 



function: 

F(VPDL) = .4 const + F°^ C (VPDL) • 4, (40) 

where afj and A cons t are the two free parameters. The 
constant asymmetry term allows for a contribution from 
possible additional background sources of asymmetry 
that have not been considered in this analysis. For this 
study we neglect any uncertainties on F^ C . The results 
are as follows: 

4 = (0.51 ± 0.86)% {jiD channel), (41) 
4 = (1.25 ±0.87)% (p-D* channel). (42) 

These values are consistent with those from the full anal- 
ysis method. The uncertainties are larger as a result of 
the additional parameter in the fit. The constant asym- 
metry parameter converges to values consistent with zero 
for both channels: 

A const = (-0.03 ± 0.23)% (\xD channel), (43) 
A const = (-0.09 ± 0.21)% (fiD* channel), (44) 

demonstrating that any possible residual background 
asymmetries not accounted for are small, as expected. 

XI. COMBINATIONS WITH OTHER 
MEASUREMENTS 

This measurement of a~, can be combined with the ex- 
isting world average from the B factories [6]. We use 
a simple weighted average, assuming that the two mea- 
surements are fully independent. The total uncertainty 



on the result presented in this article is ±0.47%, obtained 
from the addition in quadrature of statistical and system- 
atic uncertainties. The normalized weights are then 0.59 
(DO) and 0.41 (WA). We obtain: 

4 = (0.38 ±0.36)%. (45) 

This number can in turn be combined with the recent a s sl 
measurement [31], and the two-dimensional constraints 
on (flgj, a|j) from the DO measurement of the dimuon 
charge asymmetry A b si [19]. The full two-dimensional fit 
yields the following values: 

(4(comb.) = (0.07±0.27)%, (46) 
(4(comb.) = (-1.67 ±0.54)%, (47) 

where the two parameters have a correlation coefficient 
of —0.46. The results are shown in Fig. 9(a), with the 
two dimensional contours overlaid on the four constraints 
from the input measurements. The fit returns a \ 2 of 2.0 
for 2 degrees-of-freedom. The p- value of the combination 
with respect to the SM point is 0.0037, corresponding to 
an inconsistency at the 2.9 standard deviation level. 

Using only the DO measurements of a^, a^, and A^, 
we obtain the following values: 

4(comb.) = (0.10 ±0.30)%, (48) 
(4(comb.) = (-1.70 ±0.56)%, (49) 

with a correlation coefficient of —0.50. The % 2 of this fit 
is 2.9, and the standard model p-value is 0.0036, corre- 
sponding to a 2.9 standard deviation effect. Figure 9(b) 
shows the two-dimensional contours from this combina- 
tion. 



XII. CONCLUSIONS 

We have performed a measurement of the semileptonic 
mixing asymmetry from B° decays, o^, using B° — > 
jj+D^~X decays in two independent channels. We ob- 
tain <4 = [0.68 ± 0.45 (stat.) ± 0.14 (syst.)]%, which is 
consistent with the SM prediction of (-0.041 ± 0.006)%. 
The resulting precision is dominated by limited statis- 
tics in the signal channel, and is better than the cur- 
rent world-average precision obtained by combining re- 
sults from the B factories [Eq. (5)]. 

The background asymmetries are determined using 
data-driven methods, in dedicated decay channels. The 
most important background is from differences in the 
reconstruction efficiencies for positively and negatively 
charged kaons, which is of order 1%. The use of simu- 
lation is limited to measuring the relatively small (^10 
20%) fraction of signal events which do not arise from B° 
decay, and modeling the oscillation of B° mesons. 
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and NSF (USA); CEA and CNRS/IN2P3 (France); 
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FIG. 8: Final measurements of the background corrected asymmetry, in bins of VPDL(B°), for both channels. The points show 
the observed asymmetry, with the solid lines showing F^o •fflgi- Any asymmetry caused by mixing should exhibit a characteristic 
turn-on shape as the fraction of oscillated B° mesons increases. 




FIG. 9: Combination of measurements of (DO and existing world-average from B factories [6]), aj (DO [31]), and the 
two impact-parameter-binned constraints from the same-charge dimuon asymmetry A^i (DO [19]). The bands represent the 
±1 standard deviation uncertainties on each measurement. The ellipses represent the 1, 2, 3, and 4 standard deviation two- 
dimensional confidence level regions of the combination. 
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